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OVERVIEW 


Subarea Allocation Model (SAM) is designed to allocate employment, 
population, land use, and income projections for subcounty areas in the San 
Francisco Bay Region. In terms of geography, it covers the most detailed areas 
in ABAG's modelling system. Employment data includes six sectors; population 
data is presented by resident and total population, resident employment and 
number of households. SAM deals with six land use types. It projects mean 
household income and household income distribution. 


SAM is one component of ABAG's long-range projections system. The Regional 
Information System (RIS) produces regional economic, employment, demographic and 
energy projections (Brady and Yang, 1982). RIS generates employment estimates 
for 32 industry sectors. County Employment Forecasting System (CEFS), TOrecasts 
employment, income and population for counties based upon input from RIS. 
(Prastacos and Brady, 1984). The Projective Optimization Land use Information 
System (POLIS), in turn, allocates county totals to 114 superzones. Lastly, SAM 
allocates superzone statistics to 1209 census tracts. This region-county-super- 
zone-tract projection process is repeated for each projection period. 

(See Figure 1.) 


Subregional Study Areas (SRSA) and the Metropolitan Transportation 
Commission's (MTC) 550 traffic zones are two other subcounty geographic areas 
reported in SAM. SRSAs are delimited by the Local Area Formation Commission 
(LAFCO) sphere of influence boundaries for cities and aggregating census tracts 
within the sphere boundaries, including approximations of split tracts. The 
correspondence between SRSA and 1980 U.S. Census tracts is Summarized in pages 
227-232 in Projections 85 (ABAG, 1985). The correspondence between 550 zones 
and census tracts 1s available from MTC. Both SRSA and 550 zones are only 
reporting and reviewing units. 


To ensure reasonableness, accuracy and consistency, statistics produced by 
SAM are carefully reviewed at each level and in each time period. ABAG's 
Projections Technical Advisory Committee (PTAC) and county, city, and town 
planning and related agencies conduct official reviews, commenting on regional, 
county, and/or SRSA projections. All population and employment projections are 
approved by the ABAG Executive Board and become the official long-term 
projections of the agency. 


The following sections document the SAM design and its components. 
Household income projection, land use accounting model, and projection of single 
family households are discussed. In addition, a proposed procedure to improve 
the subarea projections to meet a specified development pattern is presented. 


Figure 1 ABAG PROJECTIONS SYSTEM 
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ALLOCATING EMPLOYMENT AND POPULATION 


The subarea projections model is designed to provide detailed economic and 
demographic information at the census tract level. The system is flexible and 
is designed to reflect areal variations and time fluctuations. It takes 
advantage of information on the most recent development patterns and trends on 
future growth. 


SAM utilizes ABAG's extensive land data base derived from a survey of local 
development policies (ABAG, 1986). The survey data quantified current and 
future development potentials at the census tract level, containing detailed 
information on acreage, housing units and employment. It includes residential, 
1068 | ecru and "basic" development types to define the available land 
supply. 


The allocation model projects the following six employment sectors: 


1. agriculture, forestry, and mining (SIC 01-14 excluding 074) 


2. manufacturing (SIC 19-39) 
3. wholesale ESS, RTOS Zeyh) 
4, retail (SIC 52-59) 
5. services (SIC 70-89, 074) 
6. others (balance of SIC) 


The first three sectors are considered "basic." The remaining three 
sectors are considered “local-serving" in the model. 


SAM projects employment (at place of work); households; household 
population; household income; group quaraters population; employed residents; 
land acres developed by residential, commercial, and industrial activity; and 
street and highway acreages. It has an accounting system which keeps track of 
land consumed for development during any projection period. 


The model requires the following files for its operation: 
Input: number of census tracts in each of 114 superzone files 
base year superzone file 
forecast year superzone file 
base year census tract file 
forecast year development potential file 
Output: forecast year census tract file 


The layouts of these files are given in APPENDIX A. 


The 1980 census tract file was constructed by merging (1) 1980 Census 
STF3 file for household and population data, (2) 1980 adjusted California 
Employment Development Department employment survey data, (3) ABAG 1980 
estimated land use data, and (4) 1980-2005 development potential file. 


There are 114 superzones in the ABAG region. These zones are 
aggregations of complete census tracts. Superzones were delineated to 
approximate relatively equal (in size), spatial population and employment 
subregions. The correspondence of the superzone and census tracts are 
given in APPENDIX B, 


The employment and population allocation algorithm is Specified as 
follows: 


(1) Employment 


t+1 Ps EK 


Zonal basic (local-serving) employment increment is distributed 
according to the weighted census tract share of basic (local-serving) 
employment development potentials and weighted base year employment. 


trl 

E ik -- ¢t. employment in census tract j, sector K, year t+l 
(project year) 

a -- same as above but in year t (base year) 

A Eik -- zonal change of employment between base and forecast 
year, in zone i, sector k 

Ps -- employment development potentials in census tract j; by 
basic and local-serving 

Wi, Wa -- weighting factor 


(2) Land Use 


+ ik ee bey Wie err anieMaace ree 
| 0 


Zonal basic (local-serving) land area increment is distributed 
according to the weighted census tract share of basic (local-serving) land 
development potentials and weighted base year acres. 


Gal 

Lay = et Vand acresin census’ tract 3, sector Ky year trl 
(project year) 

ae -- same as above but in year t (base year) 

Alik -- zonal change of land acres between base and forecast 
year, in zone i, sector k 

PI -- land development potentials in census tract j; by basic 
(sectors 1-3) and local-serving (sectors 4-6) 


(3) Occupied Dwelling Unit and Resident Population 


fisees ae pes ee HP; (W an ili teeta = Uae 
j j ZANHP{ AW] = pip 2 = wet 
J J 


Zonal occupied dwelling units (household population) are distributed 
according to the weighted census tract share of residential unit potentials 
and weighted base year unit. 


a 
HP ik -- c.t. occupied dwelling units (residential population) in 
census tract j, year t+l (project year) 
HP -- same as above but in year t (base year) 
Ze, -- zonal change of occupied dwelling unit (residential 
population) between base and forecast year in zone i 
pHP -- housing unit development potentials in census tract j 


(4) Group Quarters Population 


ttl ange F agp $ 
60P,. "= FAME UNE ata oae ae 
J 


Base year census tract share of zonal group quarters population 
is used to allocate zonal change of group quarters population. 


eal 
GOP 5K -- C.t. group quarters population in census tract j, year 
t+l (project year) 
Ggpt -- same as above but in year t (base year) 
ASGOP ; -- zonal change of group quarters population between base 


and forecast year in zone i 


(5) Employed Resident 


(a) find project year employed resident per dwelling unit 


Erttl 
ERt+1 ERY ‘putt! | 
So Ss =e nee * SS Sessa 
Dut+l Dus ER¢ 
ae ame 


(b) derive the projected employed residents 


toll 
ip laee 
Cy es Peo) t+] 
ERS DUS vs DU5 
(c) normalize ane so that 
ZeRtt1 = Erttl 
jd d 
ER -- employed residents 
DU -- occupied dwelling units 
1 -- zone 
J -- census tract 
t -- base year 
Seah -- projection year 


USE OF LOCAL POLICY SURVEY 


In the above allocation algorithm, the use of development potentials, in 
terms of acreage, housing units, and potential employment from ABAG's local 
policy survey data, are unique. The survey quantified the type of land, the 
timing of its availability for development, and the housing unit potentials 
(ABAG, 1986). It aggregated the specific site data to census tracts for use in 
SAM. The survey specifies supply side conditions in which household and 
employment activities are to take place. 


SAM uses development potentials as “attractors" for allocations of 
activities. At a subarea level (census tract), the greater the development 
potentials, the more the activities will be allocated there generally. 


For review of subarea allocations, development potentials are used as 
the constraint for development of acreage, housing units, and employment. 
Outputs from SAM are compared with development potentials from the local policy 
survey. 


HOUSEHOLD INCOME PROJECTION 


Projections 85 made major changes in household income projections. 
Household mean income by census tract for years 1985, 1990, 1995, 2000 and 2005 
was estimated using regression equations calibrated from the Projections 83 data 
base. Census tract mean income was then converted to ABAG's subregional study 
areas and reviewed and revised accordingly. Subregional area revisions were 
carried to the census tract data base and were reviewed again at this level of 
geography. At MTC's request, income projections were produced for each census 
tract for each of 17 income levels. 


Using a gamma function provided by A. G Wilson, mean household income was 
translated into income distribution by 17 income groups for year 2000 and 2005. 
1980 U.S. Census income distribution was compared to the gamma function 
generated distribution. Trends of 1980-2000 and 1980-2005 trans-lated gamma 
distributions were used to derive final 2000 and 2005 distributions. 
distributions. A few census tracts were hand adjusted due to their special 
characteristics. 


The following text will document the calibrated mean income projection 
regression equations, Wilson's gamma function and relevant statistics. Graphs 
of plotting on income distribution are provided. 


Mean Income Projection 


Mean household income is a variable produced by ABAG's projection project. 
The objective is to project mean household income in constant dollars by the 
smallest geographic area--1980 U.S. census tract--for every five-year period 
between 1985 and 2005. ABAG's subregional study area, MTC's 550 zone, and 
county were the other levels of geography where mean income was projected. 


Due to limitations in the available data base, only cross sectional income 
projection regression analyses were conducted in Projections 83. It was 
Strongly felt by staff that the result could be improved by incorporating time 
dimension into the analysis. Household growth over time was considered to be an 
important income predictor. Furthermore, it was felt that the results would 
improve by analyzing areas of relatively homogeneous characteristics. This was 


attempted by including county dummy variables in the Projections 83 regression 
analyses. 


By using Projections 83 census tract files for 1980 to 2000 ,2a) Set) of. 
pooled regression equations were calibrated for each of the nine Bay Area 
counties. A set of dummy variables representing area type were specified in 
the equation. The area type designation followed the results of the Social 
Area Analysis (ABAG, 1984) which classified each county in the San Francisco 
Bay Area into less than ten social areas by using factor analysis techniques. 


A large number of potential independent variables were specified and 
calibrated. They included mean household income lagged one period, growth in 
households, employed residents, land use acres by type, and land density changes 
from the local policy survey. Transformations of these variables were Specified 
too. Examples were employed residents per household in the current and previous 
period and proportions of non-residential land acres. Income growth rates for 
all periods by county, provided by Ray Brady (ABAG staff), were also specified 
in the equations. Each of the nine counties had mean income equations for years 
1980, 1985, 1990, 1995, and 2000. The base year 1980 equation was used to 
calibrate the prediction equation for other years. 


Results of the regression analyses are provided in the following tables. 
They are given by county order. Twenty-two census tracts, which are either 99 
tracts or places having predominantly institutional population, were excluded 


from the regression sampling data. Independent variables are defined as 
follows: 


ATX 


area type dummy variables in a county 


1 if the census tract belongs to area X 
QO if the census tract does not belong to area X 
X = the social area in which the regression was performed. 


ATX = 
ATX = 

where 
ERHH - employed residents per household, i.e., 


employed residents 
household 


PNR - proportion of non-residential land uses, 71.e@., 


industrial acres + commercial acres 


industrial acres + commercial acres + residential acres 


DHH - difference of household unit between current and preceding period. 
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Table 1 MEAN HOUSEHOLD INCOME REGRESSION EQUATIONS 
BY COUNTY 


Alameda County (N=284) 


Variable 1980 1985 1990 1995 2000 

AT2 4,225 -4,304 -4,446 -4,584 -4,566 
AT3 -9,783 -10,101 10,462 10,564 -10,874 
ATA =112230 ES) 125138 Si -12,644 
AT5 -13,687 -14,040 -14,693 -14,881 -15,409 
AT6 -15,081 -15,418 -16,131 -16,256 -16,789 
ATT -16,959 -17,359 -18,244 -18,428 £19,237 
ERHH 6,268 6,450 6,090 6 ,067 5 609 
PNR -4 646 -4,992 25.162 -5,153 -5,075 
DHH 0 0 0 Naa 1.75 
Constant 25,850 26,675 28,094 28,352 29,520 
R2 x 100 72.03 Bera 71.43 71.83 MP 
F 89 85 86 78 80 


Contra Costa County (N=151) 


Variable 1980 1985 1990 1995 2000 
AT2 -13,003 -13,300 -14,358 -14,460 -14,581 
AT3 ar ale 47 -21,978 -22 ,994 -23 ,503 -23 ,536 
AT4 -27,082 -27,740 -28,288 -28,983 -29,459 
AT5 -22,/05 eae aed: -23 ,880 -23,921 -24 ,354 
AT6 -30,666 = 515090 -31,809 -32,614 —3 555900 
AT7 -32,946 -33,521 -33,678 -34 404 -34,281 
ERHH 7,407 TOF 8,688 8,676 8,399 
DHH Sat 4.95 OT 6313 AS 
Constant 39,813 40,411 40,006 41,198 42,724 
R2 x 100 86.34 86.53 86.62 86.20 84.86 
F 1izZ 114 115 ket 100 


Marin County (N = 49) 


Variable 


AT2 

AT3 

AT4 

AT5 

ERHH 

PNR 

DHH 
Constant 


R2 x 100 
F 


Napa and Solano Counties 


Variable 


AT2 
AT3 
AT4 
AT5 
PNR 
DHH 
Constant 


R2 x 100 
F 


1980 
-12,418 
-16,892 
-20,349 
-16,908 

S552. 
-18,262 
-6.53 
21,803 


72.47 
15 


1980 


BoS8 
27,466 


68.27 
26 


1985 


-12,605 
-17,167 
=21,5/0 
-18,275 

i291 
-18,362 


(N = 79) 
1985 


-5,282 
=O.200 
-9,243 
-11,737 
-7,560 
coe 
20,052 


70.01 
28 
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1990 


=13,339 
-18,552 
-23,064 
-19,598 
Ly SOs 
-18,170 
-0.38 
Zo.dor 


1990 


-5,447 
ASAT) 
-9,826 
-11,246 
-7,477 


4.82 


28,488 


eae 


4) 


1995 


-5,449 
=1 997 
-10,029 
-11,634 
-7,404 
Shey 
ZO.903 


80.59 
50 


2000 


-13,833 
-19,451 
-22,808 
-21,008 
Togo 7S 
-20,316 
-4.23 
27,140 


JEN AOW 
14 


2000 


-6,081 
-8,776 
-10,284 
-12,398 
-7,506 
4.69 
30,372 


les 
41 


San Francisco (N = 145) 

Variable 1980 1985 
AT2 -14,551 -14,908 
AT3 =i Coo -17,689 
AT4 -20,/31 -21,249 
AT5 -23,586 -24,126 
AT6 -26,998 -27,616 
AT7 -27,819 -28,645 
PNR -1,827 -1,832 
DHH -5.78 -5.36 
Constant 40,910 41,841 
R@ x 100 65.77 65.59 
F 33 Sc 
San Mateo County (N - 143) 
Variable 1980 1985 
AT2 -19,906 -19,704 
AT3 -27,134 -26 ,226 
AT4 -25,120 -22,890 
AT5 -28,819 -26 ,963 
AT6 -32,466 Pa) Me hee 
AT7 -35,343 eso Lopalew As) 
ERHH Sr hlz Wet) 
DHH aoe 0 
Constant 40,243 45,116 
R@ x 100 76.18 76.48 
F 54 54 


i 


1990 1995 2000 
-15,505 -15,805 -16,058 
-18,394 -18,759 -19,061 
-22,071 -22,536 -22,873 
-25,059 -25,586 -26,037 
-28,/00 -29,571 -30,197 
-29,467 -30 ,064 -30,686 

-1,974 -1,991 -2,069 
-2.78 -2.69 -1.72 
43,414 44,407 45,185 
65.58 65.49 65.42 
oe 32 32 

1990 eds 2000 
=19,752 -20,011 -20,459 
-26,131 -26,518 -27,286 
-22,307 -22,502 -24,122 
-27,189 -27,620 -28,655 
-30,436 -30,922 -31,828 
-35,823 -36 ,076 -35,985 

6,100 6,114 alton 

0 0 20 
47,874 48,556 AT s323 
LOrou 76.47 77,06 
3) 54 50 


Santa Clara County (N = 257) 


Variable 


AT2 

AT3 

AT4 

AT5 

AT6 

AT7 

AT8 

ERHH 

PNR 
Constant 


R2 x 100 
F 


Sonoma County (N = 58) 


Variable 


AT2 

AT3 

AT4 

AT5 

AT6 

ERHH 

PNR 

DHH 
Constant 


R2 x 100 
F 


19380 


-13,446 
-13,697 
-8,955 
-15,607 
220,323 
-18,284 
-9,507 
13,399 
-7,895 
22,290 


69.83 
64 


1980 


-7,047 
-8,124 
-7,293 
-10,879 
=) C0 

7,683 
-8,263 

SH728) 

19,666 


64.49 
nh 


1985 


-13,748 
=13,053 
=919 
-15,785 
-20,417 
-18,618 
=O 59), 
13,663 
-8,372 
(ae ove) 


70.06 
64 


1985 


-6,253 
-7,637 
-6,708 
-10,095 
-5,976 
8,422 
-8,637 
6.84 
US Jee, 


[sie 
19 


1990 1995 
-14,445 -14,995 
-14,376 -14,767 
-9,570 = VOMGT 
-16,625 Agee 
21,352 Soy 
-19,466 -20 065 
-8,835 -6,744 
13,976 13,414 
Zon] -9,374 
24,038 25,995 
69.90 69.53 
64 63 

1990 1995 
-5,813 -4,818 
7890 -6 304 
-6 380 -6 885 
-10,528 -8,992 
-5,986 5,241 
8,293 7,444 
250422 -10,302 
6.09 By2il 
18,513 20,394 
67.84 60.17 
13 9 
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As expected, lagged mean household income variables were predominantly 
significant. However, these variables precluded other candidate variables 
and greatly distorted the expected signs and beta values. It was decided not 
to include them. 


Overall results are much improved from the Projections 83 income equations. 
R squares are more than doubled in most selected equations, ranging from low 0.6 
to high 0.8. Most of the dummy variables of area type are significant. Since 
most affluent areas are omitted from the equation, their signs are all negative. 
In general, areas are ranked from affluent to poor. The size of regression 
equation coefficient reflects this ordering. Specific area type effect is the 
difference between the constant and AT(i). ERHH and PNR have the expected + and 
- effects on mean income projection. Growth in households seems to have a 
negative impact on Marin and San Francisco's mean incomes. 


DHH, the difference in households (or occupied housing units) between two 
time periods is to capture the effect of household growth. In general, new 
residents are observed to have higher income potential than existing residents 
in growing areas. DHH is hypothesized to have a positive effect on mean income 
income. Jobs will bring in income; jobs in different industries and occupations 
Will produce different income streams. Since the industry/occupation of 
employed residents is not included in the projection framework, only average 
employed residents per household (ERHH) is specified. 


Many population and housing attributes of census tracts, which are related 
to income production, are not included in the ABAG projections. The area type 
classification from the social area study considered a large number of these 
attributes. The classification represents the best available characterization 
of the census tract areas. Hopefully, most of the effects of income-relevant 
attributes not specified would be captured by these area type dummy variables. 


As carried out in Projections 83, the trend of projected mean household 
income was applied to the actual base year 1980 census tract mean income to 
calibrate the forecast year mean household income model. Aggregations, of these 
census tract projected mean household incomes, by ABAG's Subregional study area 
were reviewed by Dr. Brady. Revised and rounded (to the nearest 100 dollars) 
mean household income by subregional study area was published in Projections 85. 


Census tract mean income was then normalized to the corresponding sub- 
regional study area mean household incomes. Finally, normalized tract means 
were reviewed and hand adjusted if necessary. The statistics are available in 
the census tract files. 


Mean income for MIC's 550 traffic zones was the aggregated tract mean 


weighted by the number of households. Zones are considered as a reporting 
unit. Zonal mean income was not subjected to extensive review. 
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Projection of Income Distribution 
A gamma function is chosen for the income distribution by 17 income 


groups for years 2000 and 2005 at the census tract. (See pages 137-138, A. 


G. Wilson, "Urban and Regional Model" in Geography and Planning, John 
Wiley, 1974), ; 


where 


Let 


me) 
— 
~< 
~— 
i] 


the probability that a household has income x 
T7 - the gamma function 

i - the mean income 

- the variance 

- parameters 


A value of n=5,38122 is found in ABAG's region. It is much larger 
comparing to Wilson's study (n=1.636 and 1.8). 
ex 21722/(8599)2 = 0,0002938 


n = Cx -1 = (21722)2/(8599)2@-1 = 5.38122 


In order to better understand the above distribution function, many 
plottings were made. Figure 2 plots ABAG regional income distribution with 
gamma function by three different n values. n = 5.38 curve approximates 
regional distribution in general. 


County plottings were evaluated against social Study areas. Each county 
was first plotted with 1980 census income distribution and followed by gamma 


function generated distribution. They are annotated as ‘census' and ‘gamma fn.' 
(Appendix C.) 


The last three graphs plot census distribution vs. chosen gamma function 
distribution for three census tracts in Oakland. The gamma function does not 
Simulate well in wealthy areas. This lack of fit is observed in the higher 
income ranges. Lower right corners of previous figure show this for all 
counties. (Figures 3-5.) 
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Figures 3-5 


1980 CENSUS TRACT INCOME DISTRIBUTION 
VS. CALIBRATED GAMMA DISTRIBUTION 
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Table 2 provides correlation between 1980 census and gamma function 
generated distribution for the region. Income ranges and assumed mid-point 
income for each group are also posted. 


Income distribution was projected for 1980, 2000 and 2005 by using the 
gamma function. Forecast years' (2000 and 2005) distributions were then 
adjusted to reflect the 1980 base. That is, only the trend of the gamma 
function between the base year and the forecast year is used. 


Finally, to cross check the validity of the distribution, mean income 
associated with the derived income distribution was computed for each census 
tract. This mean value was compared with the mean income obtained from the 
regression analyses. Large differences were found in tracts having relatively 
more households in the highest income group(group 17). Since group 1/7 has an 
open end range, its midpoint income was assumed to be 1.5 times of $75,000), OF 
$112,500. This may not be the best assumption. 


One plausible adjustment to the distribution was for tracts with no 
households in the 17th income group. Adjustments were made from higher to lower 
income groups. The adjustment process used an estimated aggregated income 
statistic and found appropriate number of households for assigning to lower 
income groups. This was carried out iteratively until aggregate income computed 
from the distribution and the projected mean income agreed. 


Table 2 CORRELATION OF INCOME DISTRIBUTIONS 
BETWEEN 1980 CENSUS AND GAMMA FUNCTION 


group income in 1979 ($) mid point ($) r 
1 O' =9q27500 1,667 .8699 
Z 2,900 - 4,999 4,166 .8764 
5 5,000 - 7,499 6,250 .8883 
4 7,500 -= 297999 B70 -8846 
5 10,000 - 12,499 ie 250 -8501 
6 12,500 =, 14,999 135750 S173 
7 15,000 - 17,499 16,250 ./989 
8 17 S004 195099 183750 ./884 
9 20,000 *=°22,499 DAE 250 sols 

10 22,500 - 24,999 23,750 8225 
1l 25,000 - 27,499 26,250 8342 
byl 2h 500s: 295999 25,190 sol5) 
i3 30,000 - 34,999 32,500 .8573 
14 35,000 - 39,999 37.500 .8505 
PS 40,000 - 49,999 45 ,000 Aeeiays 
16 50,000 - 74,999 62,500 7742 
17 15,000" + Fi2,500 8942 
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LAND USE ACCOUNTING 


A simple routine was developed to produce land use variables over the 
projection periods. The routine was implemented after all the variables in the 
census tract file were derived. The census tract allocation model which 
allocates all the variables from 114 Superzones to 1209 census tracts has a land 
use submodel to project tract land use. As in other submodels, the land use 
submodel takes acres by land use of the Superzones and distributes them to the 
census tracts in a given superzone. Since land use variables were not 


completely available from POLIS superzones, an independent land use accounting 
submodel was formulated. 


Land acres by the following land uses are dealt with in the procedure: 
(1) industrial developed land 

(2) commercial developed land 

(3) residential developed land 

(4) commercial and industrial available land 

(5) residential available land 


(6) street and highway developed land 


For a given census tract and in a given time period, density variables 
(employees per acre or housing units per acre) were used to transform activities 
(employment and household) to land acres. In all cases, available land statistics 
obtained from the local policy survey were the ceiling of the land development 
(j.e., change in developed land) over the projection period. 


Industrial or basic land uses include agriculture, forestry and mining, 
manufacturing and wholesale trade activities. Commercial or local-serving 
land uses include retail trade, services, and other local-serving employment 
activities. 


Initial computation of land use change was made disregarding the avail- 
ability of land. If available, local policy survey density was used. This would 
be the case if both land and "unit" statistics were available in a given tract 
in) a time? periods 9 1f local policy Survey density was unavailable, adjusted 
base year densities were used. Depending on the areal type (defined by the 
percent of developable land which were available), percents of adjustment were 
specified as follows: 
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available land 


percent available . (total land - unavailable land) 
<10% 10 - 75% TS%> 
Industrial 20 10 10 
Commercial i 3 10 
Residential 5 15 15 


Initial estimation of the developed land change was compared to the 
available land. If enough land was available to accommodate the growth, 
available lands (gross acres) were subtracted (reduced) by the growth increment, 
and developed land uses (net acres) were updated accordingly. If available land 
was less than growth, the amount of available land became the developed land use 
increase. In essence, the activities were forced into the available land by 
assuming density increases. Then the available land was set to zero for the 
next period. There was no reallocation made for either intra- or inter- 
Superzone. 


Finally, street and highway land was estimated by applying the base year 
proportion of street and highway land to the derived land increments. Indus- 
trial, commercial, and residential (developed) lands were reduced by their 
share of the street and highway land. 


In the future, if land use accounting at the superzone is working 
satisfactorily, intra-superzone reallocation of employment or occupied housing 
units due to land use availability constraints should be developed. Independent 
research on the land type and density shift by different areal types should 
be conducted to determine future land development densities. 
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PROJECTION OF PROPORTION OF SINGLE FAMILY BY 550 ZONES 


One of the projection data requirements for MIC's travel demand forecast 
is the separation of single family dwelling units from non-single (or multiple 
family) units. Single family dwelling units are defined as the sum of l-unit 
detached and l-unit attached units according to the 1980 Census definition. 
This definition differs from previous ones because it includes l-unit attached 
in addition to the l-unit detached. Townhomes and rowhouses are examples of 
l-unit attached housing units. 


Previous efforts in projection of single family housing units were 
focused primarily on residential densities. It has been observed that housing 
structural types are closely related to residential densities. Since the new 
definition of single family included l-unit attached, the calibration of the 
projection equation was expanded to include many other variables. In parti- 
cular, attention was paid to Ownership characteristics and zonal attributes. 


Based on the 1980 zonal data base, the following regression equation was 
established to project the proportion of single family units: 


P = 45.48 - 0.062 RESDEN + 6.194 ERHH - 31.28 PCTNR + 0.895E - 3 MEAN 
(R@ = 0.56, F = 101.49) 
Where P = Percent of households in single rangi yen PCL 00 
RESDEN = residential density, j.e., 


households 
residential acres 


ERHH = employed residents per household, i.e., 
employed residents 
household 
PCTNR = proportion of non-residential land, i.e. 


industrial acres + commercial acres 
industrial acres + commercial acres + residential acres 


MEAN = mean household income 


The projected percents were constrained to be non-negative and to be less 
than 100. Since the R@ was only .56, the "trend" of the forecasted percent was 
used. That is, the 1980 base year percent of Single family units is updated to 
all forecast years by multiplying the ratio of forecasted percents of year i+ 5 
Over year i. 
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PROPOSED ADJUSTMENT PROCEDURE 


Introduction 


ABAG's Projections 85 provides 1980 census tract projections in five- 
year increments. The preliminary projection, before extensive review and 
revision by ABAG staff, did not satisfactorily reflect the pattern of 
development implied by ABAG's local development policy survey (LPS). This 
necessitated time-consuming staff review and adjustments of the projection 
statistics, particularly the land use forecasts. In addition, the census tract 
allocation was made on a “yearly” basis. That is, each five-year increment 
was projected independently from other projection years. For a given census 
tract, time trends over the 25 years projection were not examined or specified 
in the model. This section will present a procedure to improve the subareal 
allocation. An experiment with the procedure and the result will be discussed. 
A general purpose adjustment routine will be provided. 


Adjustment Procedure 


RAS, as applied in adopting the national input-output coefficients as 
regional coefficients, is the type of experimental procedure used to 
improve the projection result. This procedure attempts to mechanize the 
manual adjustment process. Discussions with staff members who reviewed 
and adjusted the data acknowledged that household adjustments were carried 
out in the following fashion: 


(1) comparing growth increments of entire 25-year period (from year 
1980 to 2005) with the total residential development potential 
units from the local policy survey (LPS) data base at the 1980 
census tract; 


(2) adjusting Projections 85 according to the LPS total potential 
units. This was carried out in backward chronological order, 
from year 2005 to 1985. The time trend was smoothed except 
where local policies constrained development. In most cases, 
residential unit changes were made by trading between census 
tracts in the same subregional area. County totals were 
generally observed in the adjustment. 


(3) allowing 25-year increments to be different (in most cases 
greater) from the LPS in certain known subregional areas. These 
areas were understated in the LPS data and were permitted to grow 
beyond the extent expressed in the surveyed development 
potentials. 


Zo 


(4) adjusting population-related variables, e.g., household 
population, employed residents, and residential acres, based on 
the newly estimated residential units. Activity ratios derived 
‘prior to adjustments were used. They are household S126, 
employed residents per household, and residential density 
(households per developed residential acre). 


The proposed adjustment procedure is an iterative mechanical process which 
derives a rectangular matrix (1) to satisfy prespecified matrix row sums and 
column sums and (2) to be as similar as the distribution of the Starting base 
matrix. This procedure can potentially be applied to employment allocation 
as well as population allocation. 


As part of the ABAG projection series, there is available a detailed 
census tract projection which can be formulated as the starting base matrix. 
The projection usually incorporates recent information and revisions (as soon as 
they are available). In the top down process, there are forecasts produced by 
region, county, and subregional study areas. These statistics are formalized 
and available before the final census tract projection. They can be adapted as 
prespecified (or prescribed) column totals. Finally, total LPS derived housing 
units by census tract can be used as the row totals. For a given area to be 
Subject to the procedure, grand totals have to be identical. That is, the sum 
of column totals has to be equal to the sum of row totals. In the total 
residential unit context, total expected growth units in a projection period 
should be identical either from the county projection or from LPS. 


There are at least two points when the adjustment procedure can be applied 
at the census tract level. The first one is at the time of initial census tract 
allocation. This occurs immediately after the superzone allocation and before 
the subregional study area aggregation and review. The column totals can be 
either superzone totals or county totals. The second point is at the very final 
census tract review, i.e., before the review of census tract files by staff and 
after tract files are normalized normalized to the final Regional Planning 
Committee approved subregional study areas totals. As a matter of convenience, 
this experiment was conducted at the second point of time. 


The following experiment used (A) Projections 83 and (B) Projections 85 
increments by census tract as the starting base matrix. However, row totals 
(changes of household of 1980 to 2000) and column totals (yearly county changes 
of household) from Projections 85 were used as the targets for both experiments. 
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Experiment 


Alameda County was selected for the experiment since it represented a 
typical area with a variety of development patterns. There are about 300 census 
tracts and four periods -- 1980 to 1985, 1985 to 1990, 1990 to 1995, and 1995 to 
2000. Since year 2005 was not a projection year in Projections 83; it was not 
included in this experiment. Among the projection variables, the number of 
households is the instrumental one. Once the household units are determined, 
other population-related variables can be derived or cross checked. The 
experiment was made on the number of households. 


The ABAG LPS provides potential residential acres and units by census tract 
for the entire projection period. This information is incorporated into 
Projections 85. The differences in residential units between the year 2000 and 
year 1980 are used as the target row totals. Similarly, the four column sums 
are the desired (target) totals of the selected area. In this experiment, 
Alameda County housing unit growth increments over the four time periods were 
the target column sums. 


Result of Experiment 


Overall results were as expected. The procedure worked well in terms of 
row and column total margins convergence. This occurred at about the fourth 
iteration. This means that control totals are readily met by the proposed 
procedure. An individual cell, the household increment in a census tract over a 
five-year period, depends very much on the magnitude of the starting base 
matrix. It simulates time trend of the base matrix (either Projections 83 or 
Projections 85) closely. If reliable estimates of total growth over time in 
census tract and total growth of a county over time are available, and if a 
satisfactory starting base matrix is available, the proposed procedure will be 
very useful to generate a new growth matrix for future projection tasks at ABAG. 


One drawback of the procedure is that it cannot handle a matrix containing 
any negative cell. A negative cell means there is a decline in the activity. 
First of all, convergence is not guaranteed. In addition, if there is a decline 
in the entire projection period (a negative row sum) or if there is a one period 
decline (a negative cell) with overall increase (a positive row sum), normalized 
matrix cell values can be opposite to the intended direction of adjustment. In 
the second instance, negative cells will burden other non-negative cells of the 
same row in the procedure. This distorts the results greatly. 


Since decline of activities in certain areas is a real phenomenon, negative 
numbers have to be dealt with. As implemented in the experiment, all cells with 
negative values are set to 0 and flagged for manual adjustments. This avoids 
the convergence problem and takes care of all non-negative tracts. Since the 
procedure is intended to supplement and improve the current projection approach, 
it probably will serve this purpose well. 
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A General Purpose Routine 


An adjustment procedure subroutine is developed and is available for the 
next projection series. It is general and can be invoked to perform the 
adjustment whenever needed. Appendix D shows the annotated FORTRAN computer 
program. Figure 6 sketches the sequence of the procedure. 


Two types of applications of the procedure are anticipated; they correspond 
with two dimensions of the matrix. The first type is applicable to different 
aggregations of geography: region, county, superzone, Subregional study areas. 
The second type is applicable to different aggregations of time period data. 

It will be useful in utilizing ABAG's LPS data base. The demarcation of prior 
and post-1990 development potentials allows the procedure to be applied in two 
time periods, a useful application because short-term and long-term development 
potentials can be very different. If the procedure can treat the difference 
explicitly, the analyst will be better prepared for projection evaluation. 


An example of the procedure is provided with fictitious data. As noted in 
the computer program, the routine handles increments only. This is largely 
because the LPS provides increments of activity information. Both Starting 
and resulting matrix by increment and normal values matrix are provided in the 
example. "tb" and "td" are starting and target row sums, sums, “zb" and "zd" 
are starting and target column sums respectively. Since tract 4054 experiences 
a decline, it is signaled for special hand adjustment if needed. 
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Figure 6 FLOW CHART OF ADJUSTMENT PROCEDURE 
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APPENDIX A 


File Layout of Subarea Allocation Model 


‘can vila rs W Seas ati 
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15° =220 
Agee As) 
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33 =00 
39 - 44 
45° 2750 
5 eb 
57 - 62 
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13. =4a9 
80 - 86 
SF =293 
94 - 100 

101-7107 
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Lis ceibz? 
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Census Tract File 


Description 


County Number (Alameda =1,... Sonoma 
107 Superzone Number | 
550 MIC Zone Number 

1980 Census Tract Number 

Employment in Agriculture, Forestry, & Mi 
Employment in Manufacturing 

Employment in Wholesale Trade 
Employment in Retail Trade 

Employment in Services 

Employment in Others 

Industrial Land Acres 

Commercial Land Acres 

Street and Highway Acres 

Residential Available Acres 

Industrial and Commercial Acres 

Vacant Unavailable Acres 

Household Number 

Household Population 

Group Quarters Population 

Employed Resident 

Fraction Single Families 

Mean Household Income 

Median Household Income 


Number of Household in 1979 Income $ 
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= 9) 


ning 
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Census Tract File 


Fields Size Description 

143 - 149 wPAe Number of Household in 1979 Income $12,500 - 29,999 
150: = 156 abe) Number of Household in 1979 Income $30,000 or more 
157 - 163 F7<0 Total Acres 

164 - 170 F7.0 Residential Acres 

UW At Wane Fa a Development Potentials, Residential Acres 

178 - 184 Evra Development Potentials, Residential Unit 

Poly cor mie bl mleul Development Potentials, Basic Acre 

192 - 198 Fay sek Development Potentials, Basic Employment 

199) = 205 Seal Development Potentials, Non-basic Acre 

206 = 212 Pe Development Potentials, Non-basic Employment 


Superzone File 


Fields ce Description 
1 - 4 14 Superzone Number 
fea F7 <0 Employment in Agriculture, Forestry & Mining 
12 468 Fy <0 Employment in Manufacturing 
19-25 Pee) Employment in Wholesale Trade 
26 - 32 ree Employment in Retail Trade 
33-939 Parad Employment in Services 
40 - 46 E70 Employment in Others 
47 - 53 Pew Industrial Land Acres 
54 - 60 PuisG Commercial Land Acres 
61 = 267 Pigs) Street and Highway Acres 
68 - 74 Ee Residential Available Acres 
75 - 81 eran Industrial and Commercial Acres 
82 - 88 hen!) Vacant Unavailable Acres 
89 - 85 Fy <0 Total Population 
96 - 102 F7.U Household Population 
103 - 109 F/O Employed Resident 
110 - 116 FI .0 Household Number 
Lien ive E7.9 Total Acres 
124 - 130 r/ sO Residential Acres 


1980 Land Use File 


Fields Size Description 
Lee] 14 1980 Census Tract Number 
8 - 15 18 Total Acres 
16 - 20 [5 Residential Acres 
21 c= 25 15 Local Serving Acres 
26 - 30 IS Basic Acres 
31: = 0 [5 Street and Highway Acres 
36 - 40 15 Rebuilding/Conversion Acres 
41 - 45 iB Vacant Available Acres 
46 - 51 16 Vacant Unavailable Acres 


Partitioned Census Tract File 


Fields Size Description 
1-4 14 Subregional Study Area Number 
Dy oe 14 Census Tract Seq. # (sorted by Superzone &C.T. #) 
9 - 14 16 1980 Cenus Tract Number 
15 - 34 A20 Name of the Subregional Study Area 
35 - 42 meio’) Total Population 
43 - 50 F8.0 Household Number 
51-458 F8.0 Household Population 
59 - 66 F8.0 Employment in Agriculture, Forestry & Mining 
67 - 74 Fos) Employment in Manufacturing 
1D = 67 F8.0 Employment in Wholesale Trade 
83 - 90 F8.0 Employment in Retail Trade 
91 - 98 F8.0 Employment in Services 
99 - 106 F8.0 Employment in Others 
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APPENDIX B 


Correspondence Table of Superzone and Census Tract 
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